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TOPICS IN SPECIAL FUNCTIONS 

G. D. ANDERSON^ M. K. VAMANAMURTHY, AND M. VUORINEN 



Abstract. The authors survey recent results in special functions, par- 
ticularly the gamma function and the Gaussian hypergeometric function. 



1. Introduction 

Conformal invariants are powerful tools in the study of quasiconformal 
mappings, and many of these have expressions in terms of special functions. 
For instance, the distortion results in geometric function theory, such as the 
quasiconformal Schwarz Lemma, involve special functions. A frequent task 
is to simplify complicated inequalities, so as to clarify the dependence on im- 
portant parameters without sacrificing sharpness. For these reasons we were 
led to study, as an independent subject, various questions for special func- 
tions such as monotoneity properties and majorants/minorants in terms of 
rational functions. These new inequalities gave refined versions of some clas- 
sical distortion theorems for quasiconformal maps. The classes of functions 
that occur include complete elliptic integrals, hypergeometric functions, and 
Euler's gamma function. The main part of our research is summarized in 
|AVV5j . 

In the later development most of our research has involved applications 
to geometric properties of quasiconformal maps. However, some of the ques- 
tions concerning special functions, raised in |AVVlj . | AVV3j . and |AVV5| . 
relate to special functions which are useful in geometric function theory in 
general, not just to quasiconformal maps. In this survey our goal is to review 
the latest developments of the latter type, due to many authors [Al| [A9| . 
[AlQlj [AlQ2l lAWl IBPRll IBPR 21 iBPSl IbFI lELl iKll lldl llOl] . 

The methods used in these studies are based on classical analysis. One 
of the technical tools is the Monotone I'Hopital's Rule, stated in the next 
paragraph, which played an important role in our work |AVV4| - |AVV5| . 
The authors discovered this result in [AVV4| . unaware that it had been used 
earlier (without the name) as a technical tool in differential geometry. See 
[Chi p. 124, Lemma 3.1] or |AQVV p. 14] for relevant remarks. 



1.1. Lemma. For —oo<a<b< oo, let g and h he real-valued functions 
that are continuous on [a, h] and dijferentiable on (a, b), with h' ^ on (a, b). 
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If g' /h' is strictly increasing (resp. decreasing) on {a,b), then the functions 

g{x) - g{a) g{x) - g{h) 

h{x) - h{a) h{x) - h{b) 

are also strictly increasing {resp. decreasing) on {a,b). 

Graphing of the functions and computer experiments in general played an 
important role in our work. For instance, the software that comes with the 
book |AVV5| provides computer programs for such experiments. 

We begin this survey by discussing some recent results on the gamma 
function, including monotoneity and convexity properties and close approx- 
imations for the Euler-Mascheroni constant. Hypergeometric functions have 
a very central role in this survey. We give here a detailed proof of the so- 
called Elliott's identity for these functions following an outline suggested by 
Andrews, Askey, and Roy in |AARl p. 138]. This identity contains, as a spe- 
cial case, the classical Legendre Relation and has been studied recently in 
[KV| and [BPSV| . After this we discuss mean values, a topic related to com- 
plete elliptic integrals and their estimation, and we present several sharp 
approximations for complete elliptic integrals. We display inequalities for 
hypergeometric functions that generalize the Landen relation, and conclude 
the paper with a remark on recent work of geometric mapping properties of 
hypergeometric functions as a function of a complex argument. 

This survey does not cover recent work on the applications of special 
functions to the change of distance under quasiconformal maps. For this 
subject the interested reader may consult |AVV5j . 



2. The F and ^ functions 
Throughout this paper F will denote Euler's gamma function, defined by 



T{z) = / e-H'^-^dt, Rez > 0, 
Jo 



and then continued analytically to the finite complex plane minus the set 
of nonpositive integers. The recurrence formula T{z + 1) = zT{z) yields 
F(n + 1) = n\ for any positive integer n. We also use the fact that F(l/2) = 
^/^^. The beta function is related to the gamma function by B{a,b) = 
F(a)F(6)/F(a + b). The logarithmic derivative of the gamma function will 
be denoted, as usual, by 

The Euler-Mascheroni constant 7 is defined as (see |A2] . [TYj . [Y]) 

n ^ 

7 = lim Dn = 0.5772156649 . . . ; D„ = V - - logn. 



= 1 



Then ^-(1) = F'(l) = -7 and ^'(1/2) = -7 - 2 log 2. For a survey of the 
gamma function see [G], and for some inequalities for the gamma and psi 
functions see lAll. 
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2.1. Approximation of the Euler-Mascheroni constant. The conver- 
gence of the sequence Dn to 7 is very slow (the speed of convergence is 
studied by Alzer [A2| ). D. W. DeTemple [Dej studied a modified sequence 
which converges faster and proved 

1 1 " 1 / 1\ 

^4(;7^<^n-7<^, where - g ^ - log + - J . 

Now let 

h{n) = Rn — 7, H(n) = n'^h{n), n ^ 1. 
Since ^'(n) = —7 — 1/n + J2k=i see that 

H{n) = {Rn - -fW = (^('^) + 1/n - log(n + ^)W- 

Some computer experiments led M. Vuorinen to conjecture that H{n) in- 
creases on the interval [l,oo) from H{1) = —7 -|- 1 — log(3/2) = 0.0173. . . 
to 1/24 = 0.0416. . .. E. A. Karatsuba proved in [Kl| that for all integers 
n ^ l,H{n) < H{n + 1), by clever use of Stirling's formula and Fourier 
series. Moreover, using the relation 7 = 1 — r'(2) she obtained, for k>l, 



12fc+l 12fc+l 



-Cfc <7-l + (log/c) d{k,r)- ^ < c^, 

r=l r=l 

where 



Cfc = 7TTT7TT + 2A:'e-^ d(A:,r) = (-l) 



(12fc)! ' V > / V y 1)1(^ + 1)' 

giving exponential convergence. Some computer experiments also seem to 
indicate that ((n -|- l)/n)^i?(n) is a decreasing convex function. 

2.2. Gamma function and volumes of balls. Formulas for geometric 
objects, such as volumes of solids and arc lengths of curves, often involve 
special functions. For example, if denotes the volume of the unit ball 
= {x : I x\ < 1} in M", and if o->n—i denotes the {n — l)-dimensional 
surface area of the unit sphere S^~^ = {x:|x| = l},n^2, then 

^" " r((n/2) + 1) ' " 

It is well known that for n ^ 7 both Qn and decrease to (cf. [AVV51 
2.28]). However, neither 0„ nor Un is monotone for n on [2, 00). On the other 

hand, J^^^"'"*^"^ decreases to 6"^/^ as n ^ 00 |AVV1[ Lemma 2.40(2)]. 

Recently H. Alzer [A4] has obtained the best possible constants a, b, A, B, 
a, P such that 

n n 



n + A jn + B 



271 1^72 V 27r 

i\" 



n i2„_ii2,„+i V 



4 G. D. ANDERSON, M. K. VAMANAMURTHY, AND M. VUORINEN 



for all integers n ^ 1. He showed that a = 2/-y/7r = 1.12837. . ., b = \fe = 
1.64872... , A = 1/2, 5 = 7r/2 - 1 = 0.57079... , a = 2 - (log vr)/ log 2 = 
0.34850 . . . , /3 = 1/2. For some related results, see [KIR] , 

2.3. Monotoneity properties. In |AnQ| it is proved that the function 

(2.4) /(.) . i5^n£±il 

X log X 

is strictly increasing from (1, cxd) onto (1 — 7, 1). In particular, for x € (1, oo), 

(2.5) x^-y-i)^-^ < r(x) < x^-^ 

The proof required the following two technical lemmas, among others: 
2.6. Lemma. The function 

oo 

n — X 



^ (n + x)^ 

n=l ' 

is positive for x G [1, 4). 
2.7. Lemma. The function 

(2.8) h{x) = x2^'(l + x) -x^'(l + x) +logr(l + x) 

is positive for all x £ [1, oo). 

It was conjectured in |AnQ| that the function / in (|2.4|) is concave on 
(l,oo). 



2.9. Horst Alzer |A2j has given an elegant proof of the monotoneity of the 
function / in (12. 4p by using the Monotone I'Hopital's Rule and the convolu- 
tion theorem for Laplace transforms. In a later paper [A3] he has improved 
the estimates in (12.511 to 



(2.10) a;"{^-i)-7 < r(x) < x^(^-i)-T, X E (0, 1), 

where a = 1 - -f = 0.42278 ...,/? = i (7rV6 - 7) = 0.53385 ... are best 
possible. If X G (l,oo), he also showed that (|2.10p holds with best constants 
a = i (vrVe - 7) = 0.53385 ...,/?= 1. 

2.11. Elbert and Laforgia |EL| have shown that the function g in Lemma [2?6 
is positive for all x > — 1. They used this result to prove that the function 
h in Lemma [221 is strictly decreasing from ( — 1,0] onto [0, 00) and strictly 
increasing from [0, 00) onto [0,oo). They also showed that /" < for x > 1, 
thus proving the Anderson-Qiu conjecture |AnQ| , where / is as in (12. 4p . 

2.12. Berg and Pedersen |BPj have shown that the function / in (|2.4p is not 
only strictly increasing from (0, 00) onto (0, 1), but is even a (nonconstant) 
so-called Bernstein function. That is, / > and /' is completely monotonic, 
i.e., /' > 0, /" < 0, > 0, . . . . In particular, the function / is strictly 
increasing and strictly concave on (0, 00). 

In fact, they have proved the stronger result that 1// is a Stieltjes trans- 
form, that is, can be written in the form 

1 da it) 

7^ = c+ / — ^, x>0, 
/(x) Jq x + t 
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where the constant c is non-negative and a is a non-negative measure on 
[0, oo) satisfying 

da{t) 



< oo. 



l + t 

In particular, for 1// they have shown by using Stirhng's formula that c = 1. 
Also they have obtained da{t) = H{t)dt, where H is the continuous density 



log|r(l-t)| + (fc-l)logt ,wfc_l^^fe_l2 
(log|r(l-t)|)2 + (A:-l)27r2' ^^^^ i,fcJ,fc-i,A, 

, i = l,2,.... 



Here log denotes the usual natural logarithm. The density H{t) tends to I/7 
as t tends to 0, and a has no mass at 0. 

2.13. In "The Lost Notebook and Other Unpublished Papers" of Ramanujan 
[Ral], the Indian mathematical genius, appears the following record: 

•T(l+.) = ^(^)'{8x3 + 4.^+. + |}"^ 

where 6x is a positive proper fraction 



= 


^ = .9675 


1/12 


= .8071 i 


^7/-^2 = .3058 


2/12 


= .6160 i 


V12 = -3014 


3/12 


= .4867 i 


^9/12 = -3041 


4/12 


= .4029 I 


'10/12 = -3118 


5/12 


= .3509 i 


'11/12 = -3227 


6/12 


= .3207 i 


)i = .3359 


- 
00 — 


1." 





Of course, the values in the above table, except ^oo, are irrational and 
hence the decimals should be nonterminating as well as nonrecurring. The 
record stated above has been the subject of intense investigations and is 
reviewed in [BCK], page 48 (Question 754). This note of Ramanujan led the 
authors of tAVVSlj to make the following conjecture. 

2.14. Conjecture. Let 

G{x) = {e/xfr{l + x)/^ 

and 

H{x) = G{xf - 8x^ - 4x2 - X = ^. 
Then H is increasing from (l,oo) into (1/100, 1/30) |AVV5[ p. 476]. 



2.15. In a nice piece of work, E. A. Karatsuba |K2j has proved the above 
conjecture. She did this by representing the function H{x) as an integral for 
which she was able to find an asymptotic development. Her work also led to 
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an interesting asymptotic formula for the gamma function: 

^(x\^( , , 2 1 11 79 3539 

r(x + l) = V^ - 8x3+4x2+x + ---— + — ^ 



e7 V 30 240x 3360x2 201600x3 

9511 10051 47474887 07 ^ a mV''^ 

~ 403200x4 ~ 716800x5 + 1277337600x6 ^^^'"^^^ n+iWj > 
where A„+i(x) = O(^p^), as x — > 00, and where each is given exphcitly 
in terms of the Bernoulh numbers. 



3. Hypergeometric functions 

Given complex numbers a, 6, and c with c 7^ 0, —1, —2, . . . , the Gaussian 
hypergeometric function is the analytic continuation to the slit plane C \ 
[1, 00) of 

(3.1) F(a,6;c;z) = 2Fi(a,6;c;z)^f;^^^;^^^, |z| < 1. 

^ (c, n) nl 

Here (a, 0) = 1 for a / 0, and (a, n) is the shifted factorial function 
(o, n) = o(a + l)(a + 2) • • • (a + n - 1) 

for n = 1, 2, 3, . . .. 

The hypergeometric function w = F(a, b; c; z) in (I3.1|l has the simple 
differentiation formula 

(3.2) 4- ^(fl' b; c;z) = —F{a + l,b + l;c + 1; z). 

az c 

The behavior of the hypergeometric function near z = 1 in the three cases 
a + b < c, a + b = c, and a + b > c, a,b, c > 0, is given by 

l"> 0, C) -Lj — r(c-a)r(c-6) ' " + < c, 
S(a, &)-F(a, 6; a + 6; z) + log(l - z) 

= i?(a,5)+0((l-z)log(l-z)), 
F(a, b; c; z) = (1 - 2:)^~""^F(c - a, c - 6; c; z), c<a + b, 

where i?(a, 6) = -2j - ^{a) - ^{b), R{a) = R{a, 1 - a), i?(i) = log 16, and 
where log denotes the principal branch of the complex logarithm. The above 
asymptotic formula for the zero-balanced case a + 6 = c is due to Ramanujan 
(see [Ai], [BH] ). This formula is implied by [ASl 15.3.10]. 

The asymptotic formula (|3.3|) gives a precise description of the behavior 
of the function F{a, b;a + b; z) near the logarithmic singularity z = 1. This 
singularity can be removed by an exponential change of variables and the 
transformed function will be nearly linear. 

3.4. Theorem. |AQVV| Fora,b> 0, let k{x) = F(a, 6; a+6; l-e"^), x > 0. 
Then k is an increasing and convex function with k'{{0,oo)) = {ab/{a + b), 

r(a + 5)/(r(a)r(6))). 

3.5. Theorem. |AQVV | Given a,b > 0, and a + b > c, d = a + b — c, the 
function ^{x) = F{a, 6; c; 1 - (1 + x)"^/'^), x > 0, is increasing and convex, 
with ^'((0,00)) = {ab/{cd), r(c)r(d)/(r(a)r(6))). 
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3.6. Gauss contiguous relations and derivative formula. The six func- 
tions -F(a± 1, b; c; z), F{a, 6ib 1; c; z), F{a, b;c±l] z) are called contiguous to 
F{a,b; c; z). Gauss discovered 15 relations between F{a,b;c;z) and pairs of 
its contiguous functions |ASl 15.2.10-15.2.27], [R21 Section 33]. If we apply 
these relations to the differentiation formula ()3.2p . we obtain the following 
useful formulas. 

3.7. Theorem. For a,b,c > 0, z £ (0, 1), let u = u{z) = F{a — 1, b; c; z), 
V = v{z) = F{a, b; c; z), ui = u{l — z), vi = v{l — z). Then 

(3.8) z^^={a-l){v-u), 

dv 

(3.9) zil — z) — = (c — a)u + (a — c + bz)v, 

dz 

and 

(3.10) ^z(l - z)F{a + l,b + l;c + I] z) = {c - a)u + {a - c + bz)v. 
Furthermore, 

(3.11) z{\—z)—{uvi+uiv—vvi) = [l—a—b) \{l—z)uvi—zuiv—{l—2z)vvi)\ . 

Formulas (I3.8l) - (l3.10p in Theorem 13.71 are well known. See, for example, 
|AAR1 2.5.8]. On the other hand, formula (l3THl . which follows from (l3^ - 
M is first proved in |AQVVt 3.13 (4)]. 

Note that the formula 
dF 

(3.12) z{l - z)— = {c- b)F{a, b - l]c;z) + {b - c + az)F{a, b; c; z) 
follows from (|3.9|) if we use the symmetry property F{a, b; c; z) = F{b, a; c; z). 

3.13. Corollary. With the notation of Theorem \3.7\ if a G (0, 1), b = 1 — a < 

c, then 



UVl + UiV — VVl = u(l) 



r(c + a-l)r(c-a + l) 



4. Hypergeometric differential equation 

The function F{a, b; c; z) satisfies the hypergeometric differential equation 

(4.1) z{l - z)w" + [c - (a + 6 + l)z]w' - abw = 0. 

Kummer discovered solutions of (|4.1|) in various domains, obtaining 24 in 
all; for a complete list of his solutions see |R2l pp. 174, 175]. 

4.2. Lemma. (1) If2c = a + b+l then both F{a, b; c; z) and F{a, 6; c; 1 — z) 
satisfy (I4.1|l in the lens-shaped region {z : < \z\ < 1, < |1 — z| < 1}. 
(2) If 2c = a + b + 1 then both F{a, b; c; z^) and F{a, c; 1 — z"^) satisfy the 
differential equation 

(4.3) z{l - z^)w" + [2c - 1 - (2a + 26 + l)z^]w' - Aabzw = 

in the common part of the disk {z : \z\ < 1} and the lemniscate {z : |1 — < 

1}- 
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Proof. By Kummer (cf. [ESI pp. 174-177]), the functions F(a, 6; c; z) 
and F{a, b;a + b + l — c;l — z) are solutions of (|4.ip in {z : < \z\ < 1} and 

: < |1 — z| < 1}, respectively. But a + b + 1 — c = c under the stated 
hypotheses. The result (2) follows from result (1) by the chain rule. □ 



4.4. Lemma. The function F{a, b; c; \/l — z^) satisfies the differential equa- 
tion 

Z^{1 - Z)zw" - {Z{1 -Z) + [c-{a + b+ \)Z\Zz^}w' - abz^w = 0, 
in the subregion of the right half-plane bounded by the lemniscate 



2cos(2t?), -7r/4 < < 7r/4, z = re'^ . Here Z 
square root indicates the principal branch. 



\/l — where the 



Proof. From (4.1), the differential equation for w = F{a,b]c;t) is given 

abw = 0. 



by 

Now put t 



t(l-t)-^ + [c-(a + 6+l)t]- 



VT^^. Then 



and 



So 



dz 


t 


dt 


z dH 


1 


IE 


z 


dz 


~7' d^ ~ 




dw 


t dw 


d^w 


t^ d^w 


1 dw 


'dt ~ 


z dz ' 


dt^ 


z^ dz^ 


z^ dz 



i(l-t) 



tW 



;W 



+ 



{a + b+l)t 



w 



abw = 0. 



Multiplying through by z^ and replacing thy Z = \/l — gives the result. □ 



If wi and W2 are two solutions of a second order differential equation, then 
their Wronskian is defined to be W{wi,W2) = w\w'2 — W2w[. 

4.5. Lemma. [AARl Lemma 3.2.6] Ifwi andw2 are two linearly independent 
solutions of (14.11) . then 

W{z) = W{WI,W2){Z) ^ 



ZC{1 _ zY+b-c+l ' 



where A is a constant. 



(Note the misprint in [AARl (3.10)], where the coefficient x(l — x) is 
missing from the first term.) 

4.6. Lemma. // 2c = a + 6 + 1 then, in the notation of Theorem \3.7\ 



(4.7) 



(c — a){uvi + uiv) + (a — l)vvi 



A-z^-^X 



Proof. If 2c = a + 6 + 1 then by Lemma 142^ 1). both vi^z) and v{\ — z) 
are solutions of (|4.ip . Since W{z) = W{vi,v){z) = v'{z)vi{z) — v{z)v[{z), 
we have 

zil-z)W{z) = z{l- z){v'vi-vv[) 

= (c — a){uvi + uiv) + (2a + 6 — 2c)vvi 
= (c — a){uvi + uiv) + (a — l)vvi. 
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Next, since 2c = a + 6 + 1, Lemma 1431 shows that z'^{l — zyw{z) = A, 
and the result follows. □ 

Note that in the particular case c = l,a = b = ^ the right side of (|4.7|1 is 
constant and the result is similar to Corollary 13.131 This particular case is 
Legendre's Relation (I5.3|l . and this proof of it is due to Duren |Du| . 

4.8. Lemma. If a,b > 0,c > 1, and 2c = a + b + 1, then the constant A 
in Lemma 14.61 is given by A = (T (c))'^ / {T {a)T (b)) . In particular, if c = 1 
then Lemma 14.61 reduces to Legendre 's Relation (15. 8p for generalized elliptic 
integrals. 

Proof. The idea of the proof is to replace the possibly unbounded hyper- 
geometric functions in formula (14. Th by bounded or simpler ones. Therefore 
we consider three cases. 

Case (1): c > 2. Now a + 6>c + l>3. By (Ej]) or [ASl 15.3.3], we have 

u{z) = {l-z)'^~''F{c+l-a, c-b; c; z), ui{z) = z'^~''F{c+l-a, c-b; c; 1-z), 

v{z) = (1 — z)^~'^F{c — a,c — b;c; z), vi{z) = z^^^F[c — a,c — b;c;l — z). 
Hence 

A = {c — a)[(l — z)F{c + 1 — a, c — 6; c; z)F{c — a,c — b;c;l — z) 
+zF{c + 1 — a, c — 6; c; 1 — z)F{c — a,c — b;c] z)] 
+{a — l)F{c — a,c — b;c] z)F{c — a,c — b;c]l — z). 
Now, since a + b — c = c— 1, letting z — > 0, from (|3.3|) we get 



r(c)r(c-i) .. r(c)r(c-i) 
^ = r(a)r(5)) r(a)r(6) 

r(c)r(c-i) {T{c)f 

^ ' T{a)m r(a)r(6)' 

as claimed. 

Case (2): 1 < c < 2. Now, 1 < c < a + b < c + 1 < 3. Then 
A = {c — a)[{l — zY^^u{z)F{c—a,c—b;c; l — z) + z'^~^ui{z)F{c—a,c — b;c;z)] 

+{a — l)F(c — a,c — b;c; z)F(c — a,c — b;c;l — z). 
Now letting z — > 0, from (|3.3p . as in Case (1), we get 



(c-1) 



r(a)r(6) ^ ' r(a)r(6) 
r(c)r(c-i) {T{c)f 



r(a)r(6) r(a)r(6)' 

as claimed. 

Case (3): c = 1. Now a + b=l. Then 

A= {I - a)[u{z)vi{z) + ui{z)v{z) - v{z)vi{z)\ 
= (1 - a)ui{z)v{z) + (1 - a)vi{z)[u{z) - v{z)]. 
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From (EH Ex. 21(4), p. 71] we have 

u{z) — v{z) = F{a — l,b]c; z) — F{a,b;c; z) 



c-b 

zF{a, b]c + l;z) — zF{a, b; c; z), 



so that 

u{z) — v{z) c — b 



F{a, b;c + l;z) — F{a, b; c; z) —b/c, 



z c 

as z — > 0. Also, by (|3.3p . zvi{z) — > as z ^ 0. Hence, letting z ^ 0, we get 

r(c)r(c + l-a-6) 



A = (1 - a)ui(l) = (1 - a) 



r{c + l-a)T{c-b) 

(r(c))2 (r(c))2 



(1 - a)r(o)r(5) r(a)r(5)' 

as claimed. 

Note that, in Case (3), T{c) = T{1) = l,r(6) = r(l-a), and thus by [ASl 
6.1.17] A = l/(r(a)r(l - a)) = (sin7ra)/7r. □ 

For rational triples (a, b, c) there are numerous cases where the hyper- 
geometric function F{a,b;c;z) reduces to a simpler function (see [PBMj ). 
Other important particular cases are generalized elliptic integrals, which we 
will now discuss. For a,r £ (0, 1), the generalized elliptic integral of the first 
kind is given by 

3Ca = 3Ca(r) = ^F(a,l-a;l;r2) 

= (sin7ra) / (tant)^"^"(l - sin^ t)"" dt, 
Jo 

%', = %',{r)=%a{r'). 
We also define 

2sm(7ra) Xa{r) 
The invariant of the linear differential equation 
(4.9) w"+pw' + qw = Q, 

where p and q are functions of z, is defined to be 

(cf. [R2,p.9]). If wi and W2 are two linearly independent solutions of (|4.9|) . 
then their quotient w = W2/W1 satisfies the differential equation 

sUz) = 21, 

where S^j is the Schwarzian derivative 

yw J 2 \w' 
and the primes indicate differentiations (cf. |R2l pp. 18,19]). 
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From these considerations and the fact that Xai^) and OC'^{r) are hnearly 
independent solutions of (14. 3|) (see |AQVV[ (1.11)]), it follows that w = Ha{r) 
satisfies the differential equation 

-8a(l - a) 1 + 6r2 - 3r^ 
- + 2r2(r')4 ■ 

The generalized elliptic integral of the second kind is given by 
8,, = 8,,{r) = ^F{a- 1,1 -a;l-r^) 



= (sinvra) / (tant)^"2''(l - sin^ t) 
Jo 

£; = £;(r) = £,(/), 



2' ' 2(1 -a)' 

For a = |, 3Ca and Ea reduce to % and £, respectively, the usual elliptic 
integrals of the first and second kind, respectively. Likewise ^i/2('^) = IJ-ir), 
the modulus of the well-known Grotzsch ring in the plane |LVj . 



4.10. Corollary. The generalized elliptic integrals %a and £„ satisfy the dif- 
ferential equations 

(4.11) ^(^')'^ + (1 - 3r2)^ - 4a(l - a)r%a = 0, 

(4.12) ^(^'f^ + i^'f^ + 4(1 - afrE-a = 0, 
respectively. 

Proof. These follow from (Ol). □ 



For a = ^ these reduce to well-known differential equations |AVV5l pp. 
474-475], [BF] . 



5. Identities of Legendre and Elliott 

In geometric function theory the complete elliptic integrals X{r) and £(r) 
play an important role. These integrals may be defined, respectively, as 

Xir) = |F(i,i;l;r2), £(r) = |F( i , - 1 ; 1; r^), 

for — 1 < r < 1. These are 3Ca(r) and £a(?"), respectively, with a = ^. We 
also consider the functions 

X' = X'{r) = %{r'), < r < 1, 

3C(0) = 7r/2, D<:(l") = +oo, 

and 

£' = £'(r) = £(r'), ^ r ^ 1, 

where r' = \/l — r^. For example, these functions occur in the following 
quasiconformal counterpart of the Schwarz Lemma [LVj : 
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5.1. Theorem. For K S [l,oo), let w he a K -quasiconformal mapping of 
the unit disk D = {z : \z\ < 1} into the unit disk D' = {w : {wl < 1} with 
w{0) = 0. Then 

\w{z)\ ^ (Pk{\z\), 

where 

(5.2) ^^(^)^^-i(_^(,)) and ;.(r) = 

This result is sharp in the sense that for each z G D and K G [l,oo) there is 
an extremal K -quasiconformal mapping that takes the unit disk D onto the 
unit disk D' with w{0) = and \w{z)\ = ipKi\z\) {see |LVl p. 63]). 

It is well known jBF] that the complete elliptic integrals JC and £ satisfy 
the Legendre relation 

(5.3) e.x' + £.'%- X%' = 



For several proofs of (|5.3p see [Duj . 

In 1904, E. B. Elhott [E] (cf. fAVV3| ) obtained the following generalization 
of this result. 

5.4. Theorem. // a, 6, c ^ and < j; < 1 then 

T{a + b+l)T{b + c+l) 



(5.5) F1F2 + F3F4 - F2F3 

^ ' 12^34 Yia + b + c+l)r{b + 



where 



Fi= F{^ + a,-^-c;l + a + b;x], 
F2 = F[^-a,^ + c;l + b + c;l 
F3 = F{^ + a,^-c;l + a + b;x 
F4 = F{ -^-a,^ + c;l + b + c;l-x 



Clearly (|5.3p is a special case of (|5.5p . when a = b = c = and x = 
r^. For a discussion of generalizations of Legendre's Relation see Karatsuba 
and Vuorinen |KV| and Balasubramanian, Ponnusamy, Sunanda Naik, and 
Vuorinen |BPSV| . 

Elliott proved (|5.5p by a clever change of variables in multiple integrals. 
Another proof was suggested without details in |AARl p. 138], and here we 
provide the missing details. 

Proof of Theorem 15.41 In particular, let yi = F3, 2/2 = x~"~^(l — 
x)'>+^F2. Then by [R2l pp. 174, 175] or \KM\ (3.2.12), (3.2.13)], yi and ?/2 
are linearly independent solutions of (|4.ip . 

By dSIll), 

(5.6) xil-x)y[ = (a + 6 + c + + [-(a + & + c + ^) + (a + hxjFs, 
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and by ([3?9]l . 

x{l - x)y'2 =x(l - x) [ - (a + b)x-''-^-\l - x)''+' 



- {b + c)x-"'-\l - xy-^^-^]F2 
^^•^^ -x-''-'{l-x)'+'[{a + b + c+^)F4 

+ [-{a + b + c+^) + {c+^){l-x))F2]. 

Multiplying (15. 7|) by yi and (15. 6p by y2 and subtracting, we obtain 



So 



1, 



x(l - x)iy2y[ - yiy'2) = {a + b + c+ -)x-''-\l - xf+^F^F, 



2-- ' 
+ [ - (a + 6 + c+^) + (a + ^)x]x-"-''(l-x)^+^F2F3 

+ x(l - x)[{a + 6)x-"-''-i(l - xf^" 
+ {b + c)x'''-\l-xf+'-^]F2F^ 

+ x-''-\l - xf+''{a + b + c+ ]^)F^Fi 

+ x-"-^(l - xf^" [-{a + b + c+\) + {c+\){l-x)\ F2F- 



2' 

: (a + 6 + c + ^)x-"-^(l - x)^+'=FiF2 
+ x-"-^(l-x)''+^[-(a + 6 + c+^) 
+ (a + ^)x+(a + 6)(l-x) 

+ (6 + c)x - (a + 6 + c + ^) + (c + ^) (1 - x)] F2F3 

+ x-"-^(l - x)''+^(a + b + c+ ]^)F^Fi 

■{a + b + c + ^)x-'^-^(l - x)''+^FiF2 
+ x-"-''(l - xf+'' [-[a + h + c+D] F2F3 
+ {a + b + c+ ]^)x-''-\l - xf+^F^iFi. 



x{l-x)W{y2,yi) 



3 



a-b^l _ xf+c^^ + 5 + c + ^) [F1F2 + F3F4 - F2F3] 



by Lemma |431 Thus 

F1F2 + F3F4 — F2F3 = A, 

where ^ is a constant. 
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Now, by (lOl . 

^ ^ , ^/l 1 , X T(a + b+l)T(b + c + l) 

F,F. .end. to F(- + „, - c; „ + . + 1; l) = ^ . j^, ^ ^ ^ 3! 

as X — > 1, since ^ + a + (— ^ — c) = a — c<a + 6+l. 
Next 

-^3-^4 — F3F2 = F^{F4 — F2), 
where F4 — F2 ~ const • (1 — x)^ + 0((1 — x)'^), and 

r(a + 6 + l)r(6 + c+l) 



3 



r(6 + i)r(a + 6 + c+i; 

, ^^^_L?i_Ll 

2^2 



ifa + i + i- c<a + 5+l, or-c<6, i.e., 6 > or c > 0. If -c = 6 = 0, 



then 

1 1 



fi(a + i,i-c) 



^3 = ' ' f + 0((1 - x) log(l - x)) 



by (13. 3|) . In either case the product ^3(^4 — F2) tends to as x ^ 1. The 
third case a+^ + ^ — c > a+b+1 is impossible since we are assuming that 6, c 
are nonnegative. Thus A = T{a + b + l)T{b + c+l)/{T{b+^)T{a + b + c+^)), 
as desired. □ 

The generahzed elhptic integrals satisfy the identity 

(5-8) £.a%'a + £-'a^a - '^a^'a = -rrr^ T- 

4(1 - a) 

This follows from Elliott's formula (|5.5p and contains the classical relation 
of Legendre (|5.3|1 as a special case. 

Finally, we record the following formula of Kummer [Kuml p. 63, Form. 
30]: 

F{a, b;a + b - c + l;l - x)F{a + 1, & + 1; c + 1; x) 



c 



a+b-c+l 



F{a, b; c; x)F{a + l,6+l;a + 6- c + 2;l-x) 



z?x-(i-x)— ^=r(a + 6-c + l)r(c + l) 



r(a + i)r(6 + i) 

This formula, like Elliott's identity, may be rewritten in many different ways 
if we use the contiguous relations of Gauss. Note also the special case c = 
a + b-c+l. 



6. Mean values 

The arithmetic-geometric mean of positive numbers a, 6 > is the limit 
AGM{a,b) = lima^ = lim6„, 
where ao = a, 60 = b, and for n = 0, 1, 2, 3, 



fln+l = Mo-n, bn) = {an + bn)/2, &n+l = G(a„, 6„) = \/(hAi-, 

are the arithmetic and geometric means of an and 6„, resp. For a mean value 
M, we also consider the t-modification defined as 

Mt{a,b) = M{a\b^f'K 
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For example, the power mean of a, 6 > is 
and the logarithmic mean is 

L{a,b) = 

log(6/a) 

The power mean is the t-modification of the arithmetic mean Ai{a,b). 

The connection between mean values and elliptic integrals is provided by 
Gauss's amazing result 

This formula motivates the question of finding minorant/majorant functions 
for 3C(r) in terms of mean values. For a fixed x > the function t i— > 
Lt{l,x),t > 0, increases with t by [VVl Theorem 1.2 (1)]. The two-sided 
inequality 

L^/2{l,x)> AG{l,x)> L{l,x) 

holds; the second inequality was pointed out in |CV| . and the first one, due 
to J. and P. Borwein [BB2] , proves a sharp estimate settling a question raised 
in connection with |VVj. Combined with the identity above, this inequality 
yields a very precise inequality for X{r). 

Several inequalities between mean values have been proved recently. See, 
for instance, |XlQ2] , [QS] , [Hi], [S2], [S3|, [Tj, [C], and [E]. 

Finally, we remark that the arithmetic-geometric mean, together with 
Legendre's Relation, played a central role in a rapidly converging algorithm 
for the number vr in [Saj . See also |BBH IhI ILel ILuj . 

7. Approximation of elliptic integrals 

Efficient algorithms for the numerical evaluation of X{r) and £(r) are 
based on the arithmetic-geometric mean iteration of Gauss. This fact led to 
some close major ant /minor ant functions for X{r) in terms of mean values in 

m- 

Next, let a and b be the semiaxes of an ellipse with o > 6 and eccentricity 
e = \J — b"^ /a, and let L{a, b) denote the arc length of the ellipse. Without 
loss of generality we take a = 1. In 1742, Maclaurin (cf. |AB| ) determined 
that 

=4£(e) =27r -2^1(1, -i;l;e2). 

In 1883, Muir (cf. |AB| ) proposed that L{l,b) could be approximated by 
the expression 27r[(l -|- 6^/^)/2]^/^. Since this expression has a close resem- 
blance to the power mean values studied in [VV| . it is natural to study the 
sharpness of this approximation. Close numerical examination of the error 
in this approximation led Vuorinen |V2j to conjecture that Muir's approxi- 
mation is a lower bound for the arc length. Letting r = \/l — 6^, Vuorinen 
asked whether 
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for al l r £ [0 , 1]. 

In [BPR1| Barnard and his coauthors proved that inequality (17. 1|) is true. 
In fact, they expanded both functions into Maclaurin series and proved that 
the differences of the corresponding coefficients of the two series all have the 
same sign. 

Later, the same authors |BPR2| discovered an upper bound for £ that 
complements the lower bound in (|7.ip : 

(7.2) ^£W = 2Fi(i-i;l;r2)^(i±|:^)'^', O^r^l. 
See a lso [BPS] . 

In |BPR2| the authors have considered 13 historical approximations (by 
Kepler, Euler, Peano, Muir, Ramanujan, and others) for the arc length of an 
ellipse and determined a linear ordering among them. Their main tool was 
the following Lemma 17.31 on generalized hypergeometric functions. These 
functions are defined by the formula 

pFq[ai,a2,- ■ ■ ,ap;6i,62, • • • ,bq;z) = 1 + 2^ 775 — 77 r 

n=l ^^j=i'^'^j'"') 



where p and q are positive integers and in which no denominator parameter 
bj is permitted to be zero or a negative integer. When p = 2 and q = 1, this 
reduces to the usual Gaussian hypergeometric function F{a, b; c; z). 

7.3. Lemma. Suppose a,b > 0. Then for any e satisfying < e < 1, 

3F2(-n, o, 6; 1 + a + 6, 1 + e - n; 1) > 
for all integers n ^ 1. 

7.4. Some inequalities for OC{r). At the end of the preceding section we 
pointed out that upper and lower bounds can be found for X{r) in terms 
of mean values. Another source for the approximation of 3C(r) is based on 
the asymptotic behavior at the singularity r = 1, where 3C(r) has logarith- 
mic growth. Some of the approximations motivated by this aspect will be 
discussed next. 

Anderson, Vamanamurthy, and Vuorinen [AVV2j approximated JC{r) by 
the inverse hyperbolic tangent function arth, obtaining the inequalities 

/ \ 

, ^, 7r/arthr\ , vr arthr 

for < r < 1. Further results were proved by Laforgia and Sismondi jLS| . 
Kiihnau [Ku| and Qiu [QJ proved that, for < r < 1, 

9 ^ %{r) 



8 + r2 log(4/r')' 

Qiu and Vamanamurthy |QVa| proved that 

3Cfr) 1 

^ ' < 1 + -ir'f for < r < 1. 



log (4/r') 
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Several inequalities for 3C(r) are given in [AVV51 Theorem 3.21]. Later Alzer 
IA3I showed that 



4 log 2 ' log(4/r') 



for < r < 1. He also showed that the constants j and 7r/(41og2) — 1 in 
the above inequalities are best possible. 

For further refinements, see |QVul[ (2.24)] and |Be| . 

Alzer and Qiu |A1Q1| have written a related manuscript in which, besides 
proving many inequalities for complete elliptic integrals, they have refined 
(|7.5p by proving that 

TT / arth r \ 3/4 tj- arth r 

< X{r) < 



2 V r / ^ ' 2 r 

They also showed that 3/4 and 1 are the best exponents for (arth r)/r on 
the left and right, respectively. 

One of the interesting tools of these authors is the following lemma of 
Biernaki and Krzyz [BKj (for a detailed proof see |PVlj ): 

7.6. Lemma. Let r„ and Sn, n = 1, 2, . . . be real numbers, and let the power 
series R{x) = YlnLi ^n.x"' and S{x) = YlnLi ^n^" convergent for \x\ < 
1. If Sn > for n = 1,2,..., and ifrn/sn is strictly increasing {resp. 
decreasing) for n = 1,2,..., then the function R/S is strictly increasing 
{resp. decreasing) on (0,1). 

7.7. Generalized elliptic integrals. For the case of generalized elliptic 
integrals some inequalities are given in |AQVV| . B. C. Carlson has intro- 
duced some standard forms for elliptic integrals involving certain symmetric 
integrals. Approximations for these functions can be found in [CG| . 

8. Landen inequalities 

It is well known (cf. [BF| ) that the complete elliptic integral of the first 
kind satisfies the Landen identities 

x(|^) = ,i + rW0, :K(i^) = if:.V). 

Recall that %{r) = ^F{^, i;l;r^). It is thus natural to consider, as sug- 
gested in |AVV3| . the problem of finding an analogue of these formulas for 
the zero-balanced hypergeometric function F{a, b; c; r) for a,b,c > and 
a + b = c, at least when the parameters (o, b, c) are close to (^, ^, 1). From 
(3.3) it is clear that F{a,b; c;r'^) has a logarithmic singularity at r = 1, if 
a,b > 0, c = a + b (cf. (AAR] V Some refinements of the growth estimates 
were given in [ABRVV] and |PV1| . 

Qiu and Vuorinen |QVul| proved the following Landen-type inequalities: 
For a,b € (0, 1), c = a + b, 



F(a,b;c; ^ (1 + r)F(a, 6; c; r 



2\ 



< F{ a,b;c 



18 G. D. ANDERSON, M. K. VAMANAMURTHY, AND M. VUORINEN 



and 



F{a,b;c;l - r^) F( a,b;c; ' ' 



1 + 

F(a,&;c;l-r2) + l(i?-logl6) 



with equality in each instance if and only if a = 6 = ^. Here B = B{a,b), 
the beta function, and R = R{a, b) = —27 — ^(0) — ^(6), where ^' is as given 
in Section 2. 

9. HYPERGEOMETRIC series AS AN ANALYTIC FUNCTION 

For rational triples (a, 6, c) the hypergeometric function often can be ex- 
pressed in terms of elementary functions. Long lists with such triples con- 
taining hundreds of functions can be found in [PBM| . For example, the 
functions 

/(z) = zF(l,l;2;z) =-log(l-z) 

and 

g(z)^zF(l,l/2;3/2;z2) = ilog('i±i 

2 \1 — 2; 

have the property that they both map the unit disk into a strip domain. 
Observing that they both correspond to the case c = a + b one may ask 
(see |PVll IPV2| ) whether there exists 5 > such that zF{a, 6; o + 6; z) and 
zF{a, b;a + b; z'^) with a, 6 G (0, 6) map into a strip domain. 

Membership of hypergeometric functions in some special classes of univa- 
lent functions is studied in (PVTl [PV2l [BPV2| . 
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